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Abstract 


An analytical transformation method which can map arbitrary physical coordinate 
grid distribution into desired computational coordinate with uniform grid distri- 
bution is derived. The transformation function and its higher derivatives are 
differentiable. Salient features include; 1} precise control of grid sizes, 2) 
more than one location of clustered grids, 3) exact positioning of particular 

j? i computational nodes in the physical plane, 4) ensuring several patches of uni- 

w/ formly spaced grids in the physical plane for the higher accuracies (such as at 

the boundaries), etc., while keeping the variation of grid spacing continuous to 

avoid numerical instability. 


Work was performed while author was employed by Old Dominion University 
Research Foundation under NASA Research Grant NSG 1 087 ( J . E.' Harris, 
Technical Monitor, High-Speed Aerodynamics Division, NASA LRC) 
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TRANSFORMATION FUNCTIONS 


y = Coordinate in physical plane. 

n = Coordinate in computational plane. 
Nodes are uniformly spaced. 


Parameters Specified 

N = Number of pivot 

n (i=l » 2, — N) 
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BUILDING BLOCKS OF TRANSFORMATION FUNCTION f( n ) 


n (i = l, 2, — N) = "Pivot points," particular specifiable values of n where 
p i f'(ri) assumes local maxima. 

Y = 2.326, a convenient constant for the scaling cu . 

a. (i r l, 2, — N) = "Width parameter," specifies width in n in which 90 percent 

of grid size variation takes place around the pivot n p -j. 
a<0 specifies increasing grid sizes, and a>0 specifies 
decreasing grid size at pivot n pi -- 

e . (i=0, 1,2 — N) = step heights for the pivots, which decides the ratio between 

the sizes of node spacings in yi on both sides of the pivot. 
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THE EFFECTS OF CONTROL PARAMETERS, EXAMPLE 1 


Grids are clustered only at n m i n with a single pivot. 
Specified Parameters 

Total node points = 51 
Number of pivots = 1 

n Pl = 0.6 

a l =-0.2 

6 0 = 0.1 
B 1 = 0.4 


Both ends of physical plane nodes have constant spacing. 






EFFECTS OF CONTROL PARAMETERS, EXAMPLE 3 


Grids are clustered only at an interior point with a pair of pivots 
Specified Parameters 

Total node points = 51 
Number of pivots = 2 
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EFFECTS OF CONTROL PARAMETERS, EXAMPLE 5 


Figure demonstrates discontinuous variation of grid sizes (small 
compare to pivot spacing) with multi-pivots. 


Specified Parameters 
Total node points = 51 
Number of pivots = 2 
B 0 = 0.1 

rip ^ - 0.6, otj = -0.001, Bj 

r)n “ 0.9, do = -0.001, Bo 


ND. OF GRID = SI NO. OF PIVOT- 
5(d)- 0.000 0.300 

fl(d)= -0.001-0.001 
B(d)= 0.200 0.700 


2 B0= 0.10000 


CONSTITUTIVE RELATIONS FOR THE 
DETERMINATION OF B'S 


n c . (k=0, 1, 2, — K) = "Cluster Points," specifiable values of n which must 

coincide with particular predetermined values of 
y ck ( k=0 , 1, 2, --K) of corresponding k's. The values 
nek (hence y also) may or may not necessarily be 
node points. c k 

p, (k=0, 1, 2, — K) = Ratio between the immediately neighboring Ay minima 

and maxima at the cluster point k. 

2K simultaneous linear algebraic equations for 2K unknowns, e n *, & 9 , fo — 
are solved by the program “FIXBY". 
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APPLICATION, EXAMPLE 1 


The objectives were (1) to transform the physical coordinate y in the 

range of 0 ^ y < 2.5 to computational coor- 
dinate n in the range of 0 f n < 1 . 0 ; 

( 2 ) to obtain grids clustered at y=0 and y=l with 
the same order of minimum grid sizes and the 
same rates of increase of grid spacing Ay away 
from y=0 and y=1.0; 

(3) to use a total of 41 nodes where y= 0 , 1.0, 2.5 
which are nodes in the physical as v.'ell as in 
the computational plane (i.e., cluster points); 

(4) to have grids near y=2.5 which are fairly uniform 

Program "FIXBY" is used to compute Bg, B,, B 2 > 63 , with the specified input 
values; 

n pi = 0.2, 0.425, 0.825 
a.,- = -0.17, 0.17, -0.17 for i = 1, 2, 3 

n r = 0. , 0. 625 

y Ck = 0. , 1.0 

P k = 0.0025, 0.0015 for k = 1, 2 

The choice of np.- decides the way in which the number of nodes are 
partitioned into different regions and the choice of a-j decides the ratio 
with which the grid size varies. 

Computed B's are 

B i = 0.011071747, 4.417627, 4.4220557, 8.5569778 ( i =0 , 1, 2, 3) 


V. 
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APPLICATION, EXAMPLE 2 




The objectives were to have 


( 1 ) 

( 2 ) 

(3) 

(4) 


a total of 41 nodes ; 

variable ranges of 0 - y - 8.0 and 
0 - n - 1.0, respectively; 

one interior clustering point at y=0; and 

grids near y=0 and y=8.0 which are 
fairly uniform. 


Inputs to program FIXBY were; 
np.j = 0.3, 0.7 
af = 0.17, -0.17 (1-1,2) 
n Cl = 0.5 
ycj = 3 
pj = 0.25 
Computed B's were; 

6 i » 8.5714166, 6.4285625, 13.095257 (1=0 , 1, 2) 
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ORIGIN 




0F POOR Q’JAUT i 

NO. GRIC 

> = HI ND. 

CLU57ER= 1 ND. PIVDT= 2 

B0= B.57IHIGG39 

5(l))= 

B. 300000000 

0.700000000 


R(d) = 

0. 170000000 

-0.170000000 


B(d)= 

6.H2B5G2HB0 

13.035255990 




FLOATING GRID CLUSTER REGION 


When pivots are sufficiently separated so that there will be a region of 
"f(n) = constant" between every pivot (small |a| case) is a particular case. 

A separate program "FLOAT" is programmed for this case. Floating the location 
of clustering and varying the degree of clustering as well as locating new 
clustering regions during the computation (of a finite difference solution 
method) can be conveniently accomplished with this code. Figure shows a trans 
formation of "FLOAT" code with a particular set of input parameters. "FLOAT" 
redistributed the uniform computational grids n in the physical plane from 
y.(l) to y.(2). 
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